arXiv:hep-lat/0512041v3 27 Sep 2006 


UUITP-28/05 
ITEP-L AT/2005-29 

Embedded monopoles in quark eigenmodes 
in SU(2) Yang-Mills Theory 

M. N. Chernodub 1,2 and S.M. Morozov 1 

1 ITEP, B.Cheremushkinskaya 25, Moscow, 117259, Russia 
s Department of Theoretical Physics, Uppsala University, P.O. Box 803, S-75108, Uppsala, Sweden 

We study the embedded QCD monopoles (“quark monopoles”) using low-lying 
eigenmodes of the overlap Dirac operator in zero- and finite-temperature SU(2) 
Yang-Mills theory on the lattice. These monopoles correspond to the gauge-invariant 
hedgehogs in the quark-antiquark condensates. The monopoles were suggested to be 
agents of the chiral symmetry restoration since their cores should suppress the chiral 
condensate. We study numerically the scalar, axial and chirally invariant definitions 
of the embedded monopoles and show that the monopole densities are in fact glob¬ 
ally anti-correlated with the density of the Dirac eigenmodes. We observe, that the 
embedded monopoles corresponding to low-lying Dirac eigenvalues are dense in the 
chirally invariant (high temperature) phase and dilute in the chirally broken (low 
temperature) phase. We find that the scaling of the scalar and axial monopole den¬ 
sities towards the continuum limit is similar to the scaling of the string-like objects 
while the chirally invariant monopoles scale as membranes. The excess of gluon en¬ 
ergy at monopole positions reveals that the embedded QCD monopole possesses a 
gluonic core, which is, however, empty at the very center of the monopole. 

PACS numbers: 11.15.Ha, 11.30.Rd, 14.80.Hv 


I. INTRODUCTION 

It is generally believed [0, H] that the low-temperature (confinement) and the high- 
temperature (deconfinement) phases in QCD with realistic quark masses and vanishing 
chemical potential /i are separated by a smooth crossover which takes place at temperature 
T c ps 170 MeV. As the system goes through the crossover all thermodynamic quantities 
and their derivatives change smoothly, being non-singular functions of the temperature T. 
Therefore there is no local order parameter which can distinguish between these two phases 
at [i — 0. 

Recently it was suggested |i] that a well-defined boundary between the QCD phases at 
H = 0 can still be rigorously defined as a proliferation (percolation) transition of the so- 
called “embedded QCD monopoles” or, as we also call them, “quark monopoles”. These 
monopoles are (gauge-invariant) composite objects made of quark and gluon fields. The 
monopoles are assumed to be proliferating at infinitely long distances in the high temperature 
phase while in the low-temperature phase they are moderately dilute. Contrary to Abelian 
monopoles in compact gauge theories, the embedded QCD monopoles are in general, not 
directly associated with the confining properties of the vacuum 1 . The embedded monopoles 


1 There should be however, an indirect relation between the confining properties and the embedded 
monopole dynamics since as it is well known the confinement phenomena and the chiral symmetry are 
intimately related to each other in QCD. 
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can be considered as agents of the chiral symmetry restoration: in the low-temperature phase 
the chiral condensate should be suppressed in the cores of the embedded QCD monopoles 
while outside the monopoles the chiral condensate is suggested to be non-zero. 

The assumption that the chiral phase transition should be driven by the percolation 
of such monopoles can intuitively be understood as follows [3]. At low temperatures the 
density of the embedded monopoles is low and suppression of the chiral condensate by the 
monopole cores is negligibly small. However, as the temperature increases, the density of 
the embedded monopoles gets larger and, consequently, the chiral condensate becomes more 
suppressed. One can also look at the relation between the embedded monopole density and 
the chiral condensate from another side: with an increase of the temperature the chiral 
condensate becomes weaker, and the embedded monopoles - which are energetically unfa¬ 
vorable hedgehogs in the quark-antiquark condensates - become more populated. At some 
point the chiral condensate gets low enough for the embedded QCD monopoles to become 
sufficiently dense to start proliferating themselves. 

The quark monopole in QCD is conceptually similar to the embedded defects of the Stan¬ 
dard Electroweak (EW) model ji]. These defects are called as the Nambu monopoles 5| 
and the Z-vortices jfj The Z-vortices (if they are not closed) begin and end on the Nambu 
monopoles. In the broken low-temperature phase the value of the Higgs field is suppressed 
inside the embedded EW defects, and is asymptotically non-zero outside the defects. Ac¬ 
cording to analytical estimates [7| the zf-vortices are proliferating for long distances in the 
high-temperature symmetric phase in which they form a dense percolating network. In 
the broken phase the zCvortex network is destroyed and these objects become dilute. The 
Nambu monopoles possess similar properties |g]. Numerical simulations 0,0 show that the 
percolation transition of the zf-vortices takes place both at the region of the relatively small 
Higgs mass |2j, Mh < 72 GeV, where the phase transition of the first order, and at large 
Higgs masses jlO|, Mh > 72 GeV, where the transition is a smooth analytical crossover 
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In the condensed matter physics, an onset of percolation realized in the absence of a 
thermodynamic phase transition is usually referred to as the Kertesz transition h||. The 


Ising model in an external magnetic field provides the best known example of the Kertesz 
transition which is defined with respect to the so called Fortuin-Kasteleyn (FK) clusters 0 - 
The FK clusters are sets of lattice links connecting nearest spins in the same spin states. 
These clusters are proliferating in the high temperature (disordered) phase and they are 
short-sized and dilute in the low temperature (ordered) phase. In zero magnetic field, H = 0, 
the ordered and the disordered phases are separated by a phase transition at the Curie 
temperature T c , which coincides with the percolation transition for the FK clusters. In an 
external magnetic field the phase transition is known to be absent and the ordered and 
disordered phases are connected analytically by a crossover in the T-H plane. Nevertheless, 
the phases are still separated by a Kertesz transition line T k = Th(H) which marks the 
proliferation (percolation) transition for the FK clusters. Obviously, in the zero-field limit 
the Kertesz line meets the Curie point, lim H ->o T k (H) - T c . 

The Kertesz-type transitions often appear in the gauge theories coupled to fundamental 


mater fields. Besides mentioned cases of embedded monopoles in QCD [3( and the embedded 
defects in the Electroweak model 0,00, the Kertesz line appears, for example, in the 
compact f/(l) gauge theories 0 - The manifestation of the Kertesz line in the SU (2) Higgs 
model (which is similar to the Electroweak model) can be found as the percolation of the 
center vortices [l7f. 

The picture of percolating embedded monopoles in QCD is most probably related to 
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the percolation of the hadron clusters at high temperature and non-zero density (/i 7 ^ 0 ) 
environment, which may be realized, for example, in the heavy-ion collision experiments. In 
these extreme conditions hadrons may overlap and form clusters within which the quarks 
are no more confined. The onset of the quark-gluon plasma phase is thus associated with 
the percolation transition of the hadron clusters jl.T hfil|. 

In this paper we study basic properties of the embedded (or, quark) monopoles using 
numerical simulations in the SU(2) gauge theory without dynamical matter held. The 
monopoles are defined with the help of c-valued eigenmodes of the overlap Dirac operator. 
In Section nn we describe the structure of such monopoles in the continuum space-time. 
We show that in the SU(2) Yang-Mills theory there are three types of these monopoles 
characterized by their behavior with respect to the global axial transformations. We also 
discuss the extension of our construction to the realistic SU(3) case. In the same Section 
we provide a lattice construction of the quark monopoles, which is suitable for utilization 
in numerical simulations. In Section ED we describe results of our numerical simulations for 
the density of the embedded monopoles. In Section |TV] we discuss a relation between the 
monopole density and the spectral density of the Dirac fermions. We also discuss a relation 
of our results to the Banks-Casher formula 0 - Section m is devoted to numerical analysis 
of the structure of the chromomagnetic fields around the monopoles. Our conclusions are 
given in the last Section. 


II. QUARK MONOPOLES IN CONTINUUM AND ON THE LATTICE 

The quark monopoles in QCD are analogues of the embedded (Nambu) monopoles |3, [5] 
in the Standard Electroweak model. Here we briefly outline the definition of the embedded 
QCD monopoles following Ref. |3j. For the sake of simplicity we consider the gauge theory 
with the reduced number of colors, N c = 2. An outline of the generalization of our approach 
to the bigger number of colors is given at the end of this Section. 


A. Quark monopoles in continuum SU(2) gauge theory 

Let consider the SU(2) Yang-Mills theory with one (for simplicity) species of the fermion 
held ip which transforms in the fundamental representation of the gauge group. Using ip one 
can define the bilinear functions of the fermion held, 

& = 'tp(x)VTip(x) , T= 1 , 175 , ( 1 ) 

where t = (T 1 ,r 2 ,r 3 ) are the Pauli matrices acting in the color space and 7 ^, 75 is the 
standard set of the spinor 7 -matrices in the four-dimensional space-time. The real-valued 
composite helds £5 and £4 (with the subscripts S and A corresponding to the scalar, 1, and 
axial, * 75 , operators, respectively) are scalar and, respectively, pseudoscalar (axial) helds 
from the point of view of space-time transformations. Both these helds transform as adjoint 
three-component quantities with respect to the action of the gauge group. 

In the EW model the role of the adjoint composite held m is played by the scalar triplet 
where $ is the two-component Higgs held. The EW embedded defects can then be 
formulated in terms of the classical or asymptotic configurations of the gauge and the Higgs 
helds. To make a tight link between the embedded defects in both theories we assume from 
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the very beginning that the fermion field ^ used in the definition 0 is a c-valued function. 
It is convenient to choose the held to be an eigenmode of the Dirac operator V, 

V[A\i!; x (x) = \ijj\{x ), V[A] = y M (d M + ^r a A“) + m , (2) 

corresponding to a conhguration of the gauge fields A^(x). In our numerical analysis we 
use the massless Dirac operator with m = 0. The Dirac eigenmodes are labeled by the 
eigenvalues A of the Dirac operator. The label A will be omitted in this Section for the sake 
of simplicity. 

The axial transformations are defined by the global Abelian parameter a, 

U A ( 1) : -> e“ 75 ^, $ -> ^e iQ75 . (3) 


The color vectors is and are transforming into each other under the axial transforma¬ 
tions © as follows: 


f A \ _ / cos 2a sin 2a A / \ 

\UJ V U ) V - sin2 « cos2« )\u)' 


Using two adjoint fields (JTJ) we dehne three unit color vectors 


n s = 




ni 


js X 

I is x u\ ’ 


( 4 ) 

( 5 ) 


where (A, B) and [A x B] a = e abc A b B c are, respectively, the scalar and the vector products 
in the color space and |A| = (A,A ) 1//2 is the norm of the color vector A. The last term 
in Eq. ©, Hj, is a (normalized) vector product of the scalar and axial color vectors. The 
vector Hi is interesting because it is invariant under the axial transformations ©, ©. The 
index / in the subscript of Hi stands for ’’invariant”. 

The crucial observation is to interpret the unit vectors © as directions of the correspond¬ 
ing composite adjoint Higgs field. Thus we have three Georgi-Glashow multipletsjVip, A“), 
T = S, A, /, which can be used to construct the gauge invariant ’t Hooft tensors |18|, 

FlAnr,A) = F;M)nl- 1 -^nr(D?n T )\D?nrf, T = S,A,I , (6) 

where F“ v = + ge abc A b Al is the field strength tensor for the SU(2) gauge field 

A“, and 

H' 

(Off = 5° t d, + g<?*Al, (7) 

is the adjoint covariant derivative. The ’t Hooft tensor © is the gauge-invariant field 
strength tensor for the diagonal (with respect to the color direction nr) component of the 
gauge field, 

Al = AXr, T = S,A,I. ( 8 ) 

The current of the quark monopole of the T ^ 1 type, 

kT - = JI d Al= [ dr^-As^Hx-X^fr)), ( 9 ) 
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has a <5 like singularity at the corresponding worldline C l . The monopole worldline is pa¬ 
rameterized by the vector x M (r) and the parameter r. The quark monopoles defined 

according to Eq. © are quantized and the corresponding monopole charge is conserved (he., 
the worldlines C l are closed). 

The quark monopoles k £ carry the magnetic charges with respect to the “scalar”, “axial” 
and “chirally invariant” components of the gauge field AF, Eq. (|%|). In the corresponding 
Unitary gauges n p = 5 a3 , the quark monopoles correspond to monopoles “embedded” into 
the diagonal component (JBJ). In the gauges, where the diagonal component *4]) is regular, 
such monopoles are hedgehogs in the composite quark-antiquark fields. The corresponding 
quark condensates are characterized by the typical hedgehog behavior n p ~ x a in the local 
(transverse) vicinity of the monopoles. The fact of the existence of these monopoles in QCD 
is not a dynamical fact but rather a simple (kinematical) consequence of the existence of 
the adjoint real-valued fields o, 0. Note that there is an infinite number of equivalent 
formulations of the embedded monopoles associated with triplet isovectors which are given 
by a chiral rotation © of, say, isovector with an arbitrary angle a. The dynamics of 
these monopoles is studied below. 

Let us summarize briefly: if one has the configuration of the gauge field A M and the 
configuration of the (generally massive) quark c-field if} then the location of the embedded 
quark monopoles of all three types (“scalar”, “axial” and “invariant”) can be found with the 
help of relations ©, ©, ©, ©• The quark c-fields can be defined as a set of eigenmodes 
of the Dirac operator ©, labeled by the eigenvalue A. 


B. Quark monopoles in SU(2) gauge theory on the lattice 


The lattice construction of the quark monopoles in Euclidean QCD closely resembles a 
similar lattice construction 0 of the embedded defects in the standard model of electroweak 
interactions. Consider a configuration of the lattice gauge fields U ( x , /i) and a configuration 
of the c-valued fermion matter field T(a:). Using the fermionic field T(a;) one can construct 
the composite color fields on the Euclidean lattice, 

fs(®) = 'L t (x)r a T(a;), £X( X ) = ^(x)T a ^(x ), (10) 

which are the lattice analogues of the continuum expressions ©• Then the adjoint variables 
and Aa can be used to construct the lattice unit vectors ns, n A , and n/, in a manner 
completely similar to Eq. ©• 

The next step is to define the (un-normalized) projections of the gauge field onto the 
color directions nr = npT a : 

Vy{;x, /i) = U(x, /i) + n r (x)U(x, /i)n r (x + /t), T = S , A , / . 


The field Up behaves under the action of the gauge transformation O as a gauge held, 

Vr(x, fi) —■> Q^(a;)Vr(a:, fi)fl(x + fi). 

The lattice analogue of the’t Hooft tensor © is given H3 by the compact held 9(x, fiv) e 
(—7T, 7r] defined on the plaquette P = {x, /uz/}: 


9 r (x, /in) = arg(Tr ll + n r (a;) V r (x, fi)V r (x + fi,u)V^(x + 0, fi)V^(x,u) 


( 11 ) 
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Due to the property nr(x)Vr(x, n) = Vr(x,fi)n r(x + £l), the dehnition (fill) is independent 
of the choice of the reference point x on the plaquette P. One can show that in the Unitary 
gauge, nr(x) = r 3 , the gauge invariant plaquette function (fill) coincides with the standard 
Abelian plaquette formed out of the compact Abelian fields 9£(x,n) = arg U^ 1 (x,n). 

The singularities in the compact fields 9r correspond to the quark monopoles of the 
scalar, axial and invariant types. The quark monopoles are defined on the links of dual 
lattice *{x,(i} which are dual to the cubes c XifJj of the primary lattice: 

f(x,n) =**jr- ( 12 ) 

P&dCx,!i 


where the sum is taken over all six plaquettes P forming the faces of the cube c x>li . Equa¬ 
tion m is an analog of the standard definition |24] of the Abelian monopole in the lattice 
gauge theory of compact Abelian holds. By construction, the monopole current ( 1121 ) is 
integer-valued, j r G Z, and conserved, S*j T = 0. Here the operator 6 is the lattice diver¬ 
gence. 


C. Generalization to SU(3) gauge group 


The construction of the embedded quark monopoles can be generalized to the realistic 
case of the SU(3) Yang-Mills theory. The structure of such monopoles shares similarity with 
the non-Abelian monopoles in the 577(3) gauge theory coupled to an octet Higgs held. A 
good review of the subject on the monopole configurations in the 517(3) gauge-Higgs models 
be found in Refs, 
in more detail 21 


2 cj. Below we briehy outline the construction to be discussed elsewhere 


The generators of the 517(3) gauge group are given by eight traceless matrices T a = \ a / 2 
normalized as T rT a T b = 8 ab / 2. Here A a are the standard Gell-Mann matrices. Contrary 
to the simplest case of two colors, there are two structure constants in the 577(3) gauge 
group. The totally symmetric constant f abc is defined via the relation [T a ,T b ] = if abc T c , 
while the values of the totally symmetric structure d abc can be deduced from the equation 
{T a ,T b } = + d abc T c . These constants can also be expressed via the single relation, 

APrT a T b T c = d abc + i f abc 

Similarly to the SU ( 2 ) case (JTJ) the SU (3) model admits two composite octet vectors 


& = ^(x)TTy(x) 


r = n, iqs, 


(13) 


where we use the same notations as in the 517(2) case. This should not, however, lead to a 
confusion since the 517(3) gauge group is discussed in this Section only. 

Besides the scalar octet ££ with T = 11 and the axial octet ^ a A with T = iy 5 the structure 
517(3) gauge group allows us to build the invariant held 

Z! = f abc t b s&, ( 14 ) 

which transforms in the adjoint representation of the 517(3) gauge group. The octet £/ is 
the 517(3) generalization of the invariant 517(2) triplet held £/ = x which is used in 
Eq. (0). One can explicitly show that the composite held (1141) is invariant under the global 
axial rotations which read in terms of the 517(3) octet helds (THfl) as follows: 


/ cos 2a sin 2a \ / Cs A 
y — sin 2a cos 2a J y ^ j ’ 


C a . r? c a 

sr ““ 1 4irrGr' > 


T, R = 5, A. 


(15) 
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In addition to the asymmetric constant f abc the £77(3) gauge group possess the symmetric 
structure constant d abc . Therefore one can define three additional octet fields, £g S , £ AS = 
and £ AA which form a symmetric rank -2 tensor field, 

^r 2 = ^ c ^ 2 , r 1)2 = S,A, (16) 

with respect to the global axial rotations (USD: 


Ar 2 


-Rrir^-R^r^nr' 


(17) 


Summarizing, in the realistic case of three colors we have six 2 independent structures 
which are classified with respect to the global axial rotations inn as the scalar (£“), vector 
(£g and £^), and rank -2 symmetric tensor (£g s , £° SA , and Caa)- All these structures behave 
as octet helds with respect to the £7/(3) gauge transformations. 

Note, that from a kinematical point of view any of the six octet helds can equivalently be 
used to construct the embedded quark monopoles in the £77 (3) gauge theory. The relevance 
of one or another octet held is to be figured out with the help of the dynamical considerations 
(he., with the help of numerical simulations). In the rest of this Section we describe the 
kinematical construction of the embedded monopoles in the theory with the £77(3) gauge 
group and therefore use a generic notation £ a for any of the six composite helds, or for a 
linear combinations of them. 

In order to characterize the embedded monopoles in QCD it is convenient to work in the 
Cartan-Weyl basis [22], The Cartan subgroup of the £77(3) gauge group is generated by 
two diagonal generators 

A ° °\ i / 1 0 °\ 

Hi = T 3 = - 0 -1 0 , H 2 = T 8 = ——= j 0 1 0 , (18) 

z \0 0 0/ 2 v3 ^0 0 -2 J 

which form the two-component vector H = (Hi,H 2 ). In any point of the space-time the 

octet held can be gauge-rotated to the Cartan subalgebra 


$(x) -»■ |$(x)| (h(x),H), h 2 = l. 


(19) 


where (a, b ) denotes the scalar product in the Cartan space, and h is the unit two-component 
vector pointing to the direction of the composite octet held £ in the Cartan space. 


The magnetic charge of the £77(3) monopole 0 


9m = — M* + n 2 P 2 
9 v 


( 20 ) 


is given in terms of the dual simple roots {3* with a = 1,2. The dual roots /3* are expressed 
in terms of the original simple roots (3 a of the £77(3) group as /3* = /3 a /\/3 a \ 2 . These roots 
are often chosen in the self-dual form, 


A = (1,0), P 2 — (—1/2, V3/2). 


( 21 ) 


2 In a multiflavor case there are six octet structures per one flavor, as well as a number of the flavor-mixing 
structures composed of quarks of different flavors. 





The generalization of the standard Dirac quantization of the monopole charge to the case 
of the SU (3) monopole reads as 


e ig(9M,H) _ j 


( 22 ) 


The Dirac condition ED is satisfied by Eq. ED provided the numbers n\ and n 2 are integer. 

The classification of the 517(3) monopoles qualitatively depends on the direction of the 
local Higgs field in the Cartan subalgebra. If the vector h in Eq. PI) is not orthogonal to 
any of the simple roots f3\ and /3 2 , Eq. ED, then the pattern of the symmetry breaking is 
“maximal’' 


imp! 

n 


SU( 3) -»• U(l) x U(l), 


(23) 


and the corresponding vacuum manifold is characterized by a non-trivial second homotopy 
group, 


( SU{ 3) \ 
V^i7(l) x 17(1)^ 


Z 2 . 


(24) 


Thus, the monopoles are described by the two integer numbers. These numbers are ri\ and 
n 2 which enter the definition of the monopole charge ED- 

If the asymptotic Higgs field is orthogonal either to the simple root /3\ or to the simple 
root (3 2 , then the symmetry breaking is “minimal” 


SU{ 3) -»• U{2 ). 


(25) 


This pattern also possess a non-trivial second homotopy group, 


vr 2 


( sum 

V U( 2) ) 


Z, 


(26) 


and the monopoles are characterized by one integer number which is either ri\ (if (/? 2 , h) = 0) 
or n 2 (if (/3i, h) = 0). 

The position of the SU(3) monopole can locally be determined with the help of the 
’t Hooft tensors similarly to the SU( 2) case El)- In the Cartan space the monopole is 
a point-like source of the magnetic field which is described locally as B ~ (gM, H)r /r 3 
provided r is close to the center of the monopole. The type of the symmetry breaking 
pattern corresponds to the type of the embedding of the 577(2) HP monopole into the larger 
S77(3) group. In the general case of the SU(N) gauge group the embedded monopoles can 
be formulated similarly to the case of the SU(N) gauge-Higgs monopoles reviewed in detail 
in Ref. 20 . 

It is known j2Qj that in the case of the monopole-like solutions to the classical equations 
of motion of the 577(3) gauge-Higgs model, the choice of either maximal (12311 or minimal 
ED symmetry breaking patterns depends on the details of the Higgs potential [20]. In the 
case of QCD one can realize both patterns simultaneously. Let us take any of the six octet 
fields ED , ED or ED, an d project it onto the Cartan subgroup . The color direction of this 
field in any point is usually not perpendicular to any of the root vectors ED apart from rare 
degenerate cases. Thus, if we choose any of the above octets as the effective Higgs field £, 
then it is likely that the maximal embedding pattern ED is realized. On the other hand, a 
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linear combination of any of the two octets can always be fine-tuned (again, apart from rare 
degenerate cases) to be perpendicular to a chosen root vector, so that the minimal embedding 
pattern (ESI) can also be realized in QCD. A discussion on the kinematical construction and 
well as on the possible dynamical significance of such monopoles in QCD will be published 
elsewhere jllj]. Below we discuss numerical signatures of the embedded quark monopoles in 
the simpler case with two colors. 


III. DENSITY OF QUARK MONOPOLES AT ZERO AND FINITE 

TEMPERATURE 

In order to study basic properties of the embedded QCD monopoles we perform a simula¬ 
tion of the pure SU(2) Yang-Mills model on the lattice at zero and finite temperatures. The 
technical details of numerical simulations are given in Appendix El and below we discuss 
the results of the simulations. 


A. Monopole densities and the effect of temperature 

In Figures [U (a), (b) and (c) we show the lattice densities pi att of, respectively, scalar, 
axial and invariant embedded monopoles at zero temperature. The densities, plotted in 
the units of the lattice spacing a, are shown as functions of the Dirac eigenmode energy 
A for three values of the coupling (3. Apart from a few irregular points (which we ascribe 
to statistical fluctuations), the densities are smooth functions of the eigenmode energy A. 
Moreover the scalar and axial densities at zero-temperature are decreasing functions of A for 
A > 150 MeV. The density of the chirally invariant quark monopole is a decreasing function 
for all considered values of A. The density of the chirally invariant monopoles is higher than 
the scalar and axial monopole densities. 

According to Figures E (a) and (b) the scalar and axial densities coincide with each 
other within the error bars. This observation does not contradict the statement that the 
chiral symmetry is broken in the low-temperature phase. In order to illustrate this fact 
let us consider, as a toy example, a lattice model describing the SO (2) global scalar held 
0 = (0i,02) t with two real-valued components 0i and 02- The global SO(2) symmetry 
treats the held 0 as a two-component vector in an internal space. This symmetry is an 
analogue of the axial symmetry ©, while the components 0i and 02 can be associated with 
the scalar, and axial, £4, triplets, respectively. 

Apart from a kinetic term, (c00) 2 , a generic Lagrangian of the toy model should contain 
the 5'0(2)-invariant potential 

U(0 2 ) = p 2 0 2 + A(0 2 ) 2 , (27) 

with 0 2 = 0 2 + 02 - I 11 the broken phase (/ 1 2 < 0) the scalar held develops a non-zero 
expectation value in the infinite volume system, ( 0 ) 7 ^ 0 . As a consequence, the vacuum of 
the model becomes non-invariant under the SO(2) transformations. In the symmetric phase 
(p 2 ^ 0 ) the global symmetry is restored since ( 0 ) = 0 . 

However, in a finite volume the expectation values of the components of the Higgs helds 
should vanish in both phases, (0i)y = 0 and (0 2 ) 4 / = 0 (here the subscript V in (.. ,)y 
indicates that the averages are taken in a finite volume). This happens because in the 
averages the finite number of integrals over the scalar helds 0 i ,2 (^) includes also integrations 
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FIG. 1: (Color online) The densities of the (a) scalar, (b) axial, and (c) chirally invariant embedded 
monopoles in confinement phase for (3 = 2.3493, 2.418, and 2.5 on vs. the Dirac eigenvalue A. (d) 
The same but for the deconfinement phase at T = 1.15T C . The densities are given in the units of 
the lattice spacing a. 


over all orientations of the fields in the internal space. Thus, in a finite-element system all 
the internal directions are formally equivalent. Moreover, one gets {<fri)v = {<Pi)v both in 
broken and unbroken phases because of the same reason. This is precisely the same what 
we could expect for the expectation values of the triplet fields, (£f)y = <&> v, in a finite 
volume system regardless if the chiral symmetry is broken in the thermodynamic limit or 
not. 

As a consequence, finite volume expectation values of any operator 0(£s) should be 
equivalent to the expectation value of its axial counterpart 0 (£a)- In other words, the 
choice of the axial isovector £4 in a role of the adjoint composite Higgs field is as good as the 
choice of the scalar isovector Thus the densities of the scalar and axial quark monopoles 
should coincide with each other in the finite volume. 

Apart from the quenched case, in the real QCD with dynamical quarks the breaking of 
the chiral symmetry must explicitly be seen in the densities of the embedded monopoles: 
the density of the scalar and axial monopoles must in general be different. For example, it is 
expected [3:1] that at sufficiently high temperatures the density of the axial quark monopoles 
should be higher than the density of the scalar monopoles. We expect that the effect should 
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also be observable in the finite volume because of an explicit breaking of the axial symmetry 
by the fermion determinant. In fact, the measure of the integration over the fermion fields 
is axially anomalous [ 23 }], and therefore the expectation values of the operator 0(£s) should 
be different from its axial analogue 0(£a)- In terms of the toy 50(2) model the effect of 
non-invariant measure can be emulated by an addition of an extra 50(2) breaking term (say 
SV oc ) to the potential (EH). 

In Figure |U (d) we show the density of the quark monopoles at the temperature T = 
1.15 T c corresponding to the deconfinement phase. Similarly to the zero temperature case, 
the density of the scalar and axial monopoles coincide with each other. The invariant 
monopoles are denser than the scalar/axial monopoles for all values of eigenvalue A. The 
monopole density is independent of the eigenvalue in the region 0 ^ A < 400 MeV. In the 
limit A —> 0 the densities of the scalar/axial and invariant monopoles in physical units are, 
respectively, 


lim pp hys (A) 

A—>0 r 1 


~ (3 fm) 3 [scalar and axial] 

ps (4fm) [invariant] 


T = 1.15 T c . (28) 


In the region A > 400 MeV the density of the monopoles of all three types quickly drops 
down. This observation will be confronted with the fermion spectral function in Section If VI 
To estimate the effect of temperature on the monopole density it worth comparing the 
lattice monopole densities at zero temperature for (3 = 2.3493 (shown by filled circles in 
Figures [U a),(b), and (c)) and at T = 1.15 T c for (3 = 2.35 (shown in Figure d[d)). The 
selected values of the lattice coupling [3 are very close to each other and therefore they 
correspond to almost the same value of the lattice spacing a according to Table |T] of the 
appendix El In a wide region of the Dirac eigenvalues, 0 < A < 500 MeV, the density of 
the scalar and axial monopoles at T = 0 case are noticeably smaller then the density of the 
these monopoles at T = 1.15 T c : 


Ps‘a( t = 1-167;) 

Ps“a(T = 0) 

The effect of temperature on the invariant quark monopoles is milder compared to the 
scalar/axial monopoles: 



pf x {T = 1.15 T c ) 
p\ att (T = 0) 


~ 1.5.. .2. 


(30) 


The difference in ratios ESI) and (13U1) can probably be explained by the fact that the 
chirally invariant embedded monopoles are, by definition, explicitly invariant under the ax¬ 
ial transformations ©, while the scalar and axial monopoles are not. Thus the invariant 
monopole is less insensitive to the effects of the chiral symmetry breaking/restoration com¬ 
pared to the scalar and axial monopoles. 

Summarizing, the results of this Section show that the density of the quark monopoles is 
an increasing function of the temperature in agreement with general expectations |3']. 


B. Scaling towards continuum limit 

An extrapolation to the continuum limit of numerically calculated quantities is one of 
the most important issues of the lattice simulations. In general, the monopole densities can 
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be extrapolated to the continuum with the help of the following polynomial formula: 

p latt (a) = C + v ■ a + s ■ a 2 + p ■ a 3 , (31) 

where C, v, s and p are the fitting coefficients. The terms of the order 0(a 4 ) and possible 
logarithmic corrections are neglected in Eq. m ■ Naively, if the monopoles are physical 
objects which form a gas-like ensemble then one could expect that the coefficient p - repre¬ 
senting the physical density of the monopoles - is to be non-zero while the other coefficients 
in Eq. m are vanishing. Below we show numerically that this is not the case. 

We found that the scaling of densities for all non-zero modes, A 7 ^ 0, is universal in a 
sense that the form of the scaling function does not depend on A while being sensitive to 
the monopole type. For illustrative properties we take here the eigenvalue A = 235 MeV. 
We show in Figures 0 (a) and (b) the densities of, respectively, the scalar and invariant 
embedded monopoles vs. the lattice spacing a. Since the scalar and axial monopoles have 
the same (within error bars) densities we show the data for the scalar monopoles only. In 
the same figures we show the best fitting curves for the (truncated) fitting function (Effll . 




FIG. 2: (Color online) The extrapolation of the scaling coefficients for the densities of the (a) scalar 
and axial, and (b) chirally invariant quark monopoles using various fits for A = 235 MeV. 

As it seen from Figures El the expected fit pi att oc a 3 does not work for all types of 
monopoles. The corresponding quality of the fit is \ 2 /d.o.f . = 20 and 82 for scalar/axial and 
invariant monopoles, respectively. However, the fits p latt (a) = C+v-a and p latt (a) = C+s-a 2 
give reasonable values for y 2 /d.o.f. (of the order of unity), while the coefficient C is consistent 
with zero within error bars in all our fits. Setting C = 0 we obtain that the best fits for the 
scalar/axial and invariant monopole densities is achieved, respectively, by the functions 

= s s,a(A) ■ a 2 , p l ? u (a, A) = v f (A) ■ a, A 7 ^ 0 . (32) 

In all these cases y 2 /d.o.f. ~ 1. Note that the density of the scalar/axial and invariant 
monopoles can not be well fitted by the linear and, respectively, quadratic functions of a 
since in these cases the quality of fits is as large as 10 — 20. All discussed fits are shown in 
Figures El (a) and (b) by lines. 

The similar analysis can be performed for the zero mode, Figures QKa) and (b). We find 
that the best fit functions are 


= s s,a(0) • a 2 , p\ att (a) = (7/(0) + s/(0) • a 2 , 


A = 0, 


(33) 
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FIG. 3: (Color online) The same as in Figure (J2J) but for the zero mode, A = 0. 


with x 2 /d.o.f. ~ 1.3 and 0.5, respectively. The best fit parameters are 

s SiA (0) = 1.66(3) fm 2 , (7/(0) = 0.094(2), s 7 (0) = 4.6(2) fm 2 . (34) 

Thus, the scaling properties of the density of the invariant monopoles at A 7 ^ 0, Eq. (1321) . 
and at A = 0, Eq. (El, are different from each other even on the qualitative level. 

We also attempted to determine the scaling behavior of the densities using the power 
fit of the form C a“, where C and a are fitting parameters. For the non-zero modes we 
typically get as,A ~ 2 and ay ~ 1 which is in agreement with the best fit function used 
above. For example, for the case of A = 235 MeV we get as ,a = 1-9(1) an d a/ = 1.04(7). The 
scalar/axial monopoles constructed from the of the zero Dirac mode give 05 ^( 0 ) = 2.1(1). 
The fit by the same dependence of the invariant monopole with A = 0 give a/(0) = 0.78(6) 
with higher values of x 2 /d.o.f. ss 3. Therefore the constant term Cj 7 ^ 0 in the corresponding 
fitting function (the middle formula in Eq. (1331) 1 is essential. 

The scaling coefficients ss,a and vj obtained with the help of extrapolation El to the 
continuum limit are shown in Figures 01 (a) and (b), respectively, as functions of the Dirac 
eigenvalue A. The scaling coefficient ss,a( A) of the scalar and axial monopoles has a peak 
around A ~ 150 MeV while the scaling coefficient u/(A) of the invariant mode is a monoton- 
ically decreasing function for all studied eigenvalues A. 

The behavior of the scaling coefficients s$ : a and Vj has some particularities. For example, 
we find that these scaling coefficients can be described by the formulae 

4A(V = ‘1,a ■ (ni- 7 ) > for A > 250MeV, (35) 

vf*(A) = Bi(l — —^ , for all A , (36) 

where the data for s$,a is compared with the fitting function El only in the region of large 
A with A > 250 MeV since in the small A-region the behavior of this quantity is statistically 
unclear (there is however, a noticeable tend of s$,a to decrease as A —> 0). We get (with 
X 2 /d.o.f. ss 2) the similar values for scalar and axial monopoles: 


75 = 0.54(4), — 0.15(2) fm, 

7 a = 0.44(4), l A = 0.21(2) fm. 


(37) 

(38) 








14 



FIG. 4: (Color online) The scaling coefficients for the densities of the (a) scalar and axial, and (b) 
chirally invariant quark monopoles. The fits by functions (135j) and flUD are shown by solid and 
dashed lines. 


This result suggest that the scaling exponent 7 may be close to 1/2 for scalar and axial 
types of the quark monopoles. Setting 7 s,a = 1/2 one gets 


l s = 0.171(1) fm 
l A = 0.174(1) fm 


for 7 s ,a = 1/2 • 


(39) 


The last fits are shown in Figure EI(a) by the solid and dashed lines for the scalar and axial 
monopoles, respectively. 

The scaling coefficient ty is compared to the fitting function (EHD in the whole available 
region of the eigenvalues A. We get the following best fit parameters: 

Bj = 1.09(1) fm 1 , A/ = 1.53(4) GeV. (40) 


The corresponding fit is shown in Figure 01(b) by the solid line. One can see that the scaling 
of the coefficient ty for the invariant monopoles towards small values of the eigenvalues A 
is a smooth linear function over the whole region of studied eigenvalues A. The A —> 0 
limit for the coefficients Ss,a corresponding to the scalar and axial monopoles are known less 
accurately, as it can be seen from Figure 0{b). Summarizing, in the A —> 0 limit we find: 

lim ss,a{ A) ~ 1.9(3) fm -1 , lim 17 (A) = f?/ = 1.09(1) fm -1 . (41) 

A—^0 A—>0 

As it is seen from Eqs. fSJ), (ED the scaling coefficients ss,a at A = 0 for scalar and axial 
modes seems to coincide with the corresponding limits, liniA^o ss,u(A) ~ ss,u(0). On the 
other hand, the scaling coefficient vj for the invariant monopole has a discontinuity at A = 0, 
iy(0) 7 ^ lim A ^ 0 ty(A), since the corresponding scaling formulae, Eqs. (1321) . (J23D, are different 
from each other even on the qualitative level. 


C. Cluster structure of the monopole ensembles 

The ensembles of the trajectories of the embedded monopoles can be characterized by 
percolation properties. As it happens in the case of the FK clusters in the Ising model, a 
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general ensemble of the monopole trajectories consists of clusters of different types. If in 
the thermodynamic limit at certain physical conditions there exists a non-zero probability 
to find a cluster of infinite length, then the objects are said to be percolating and are often 
called as “condensed”. In the finite volume the role of the percolating cluster is played by 
a monopole cluster with the size of the order of the system volume. Using the standard 
terminology we call the percolating clusters as “infrared” (IR) and the short-length clusters 
are referred to as “ultraviolet” (UV). 

In our studies we have used the following definitions of the IR and UV clusters '25]: 

• The largest cluster is called the IR cluster; 

• The wrapped cluster is also called the IR cluster. More precisely, for each monopole 

cluster C we calculate the sum 5' /t = ^ this sum is nonzero then the cluster 

is called the IR cluster; 

• Other clusters are called the UV clusters. 

As an example, we show in Figures 03] a), (b) and (c) the total monopole density along 
with the density of the quark monopoles in the IR and the UV clusters for scalar, axial 
and invariant monopoles, respectively. The densities are shown at zero temperature (for 
(3 = 2.3493) as functions of the eigenvalue A. One finds that the most part of monopoles 



AfMeV] A[MeV] A[MeV] 

(a) (b) (c) 

FIG. 5: (Color online) The total, infrared and ultraviolet densities of the (a) scalar, (b) axial, and 


(c) chirally invariant embedded monopoles vs. the Dirac eigenvalue A. 

of all types belongs to the IR clusters. The IR monopole density is about 3/4 of the total 
monopole density in the case of the scalar and axial monopoles, while in the axially invariant 
case almost all (about 95%) monopoles are residing in the IR clusters. Another interesting 
feature of the monopole density spectrum is that the UV part of the monopole clusters is 
almost insensitive to the value of the Dirac eigenvalue A. The UV density of the invariant 
monopoles is very small being slightly increasing function of the eigenvalue A. 

In order to estimate how the temperature affects the monopole densities, we show in 
Figure ® the densities of the monopoles in the IR and UV clusters in the deconfinement 
phase at T — 1.15T C . The lattice spacing for the data shown in Figures 0 and El is chosen 
to be almost the same. One can clearly see that the basic features of the cluster structure 
in the deconfinement phase are similar to those in the confinement phase except for the 
quantitative difference: in the deconfinement phase a bigger (compared to the confinement 
phase) fraction of the monopoles belong to the IR cluster. 

The scaling of the individual contributions (total, IR and UV) towards continuum limit 
is especially interesting. We found that the total and the IR parts of the lattice density of 
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A,[MeV] A[MeV] A[MeV] 

(a) (b) (c) 

FIG. 6: (Color online) The same as in Figure Q but for the deconfinement phase at T = 1.15T C . 


the chirally invariant monopole scale towards the continuum limit proportionally the lattice 
spacing a for all non-zero (A > 0) modes. The UV part of the lattice density does not depend 
on the coupling a at all, which indicates that this part is a lattice artifact. Note that the 
last two observations do not contradict to each other in the sense of the numerical fitting 
since the constant UV part is very small and is almost consistent with zero. For example, 
at A = 235MeV we have p l ffj v = 0.006(5). Therefore the scaling of the IR and the total 
parts should numerically be the indistinguishable from each other, and both should follow 
the functional dependence (1321) . 

As for the scalar and axial monopoles, their IR and UV clusters also scale differently. 
The IR part of the scalar and axial monopole densities (written in the lattice units) scales 
as a 2 , similarly to the total density (E2D- As for the UV part, we have found that the 
scaling of the corresponding lattice density is proportional to a. This is drastically different 
from scaling of the total and infrared parts. Unlike the invariant monopole case, in the 
case of scalar/axial monopoles there is a substantial part of the monopoles residing in the 
UV clusters. Therefore the different scaling of the UV part cannot in general be neglected. 
Unfortunately, the accuracy of our data is such that the truncated fit m with two fitting 
parameters v, s and with C = p = 0 can not give a reliable estimate of the coefficient v. In 
order to get this coefficient with a good accuracy, we fit the data for the monopole density 
in the ultraviolet clusters using the linear formula 

Pr^uv = v r V(1 > T — S, A . (42) 

An example of this fit is shown in Figure Ufa) and the corresponding coefficient of propor¬ 
tionality u uv is plotted in Figure Qb). 

Summarizing, the scaling laws for the total, IR and UV monopole densities corresponding 
to non-zero Dirac eigenvalues, are 


latt, total/ \\ 

Ps,A ( a ’ A 

= s%X(\)-a 2 , 

latt, total/ \\ 

Pi (a, A) 

= Vj otal (A) ■ a , 



Ps,A ( a > A) 

= sgA A) -« 2 , 

Pi (ffi A) 

= v} R (A) • a , 

A 7 ^ 0. 

(43) 

latt,UV/ \\ 

Ps,A ( a > V 

= (A) • ^ , 

Pi (ffi A) 

= CT(A).a°, 




As for the zero mode, the UV part of the density of the invariant monopoles is consistent 
with zero while the IR part coincide with the total monopole density within error bars. 
In the case of the scalar and axial monopoles both total, IR and UV parts satisfy the 
quadratic scaling law Moreover, Eq. (H31) indicates that in the continuum limit the 

most part of the scalar, axial, and invariant monopoles corresponding to non-zero Dirac 
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FIG. 7: (Color online) (a) Example of the extrapolation of the UV fraction of the scalar and axial 
monopole densities at the eigenvalue A = 235 MeV. The fit is done by the linear formula (1421) . (b) 
The scaling coefficient v uv vs. A for the UV part of the scalar and axial monopole densities. The 
quantity v u1 is extrapolated the continuum limit. 


eigenvalues (A ^ 0) resides predominantly in UV monopole clusters. This is not the case for 
the exact zero mode (A = 0) which even in the continuum limit may possess both IR and 
UV components of the densities. So, the exact zero modes and the non-zero modes have, in 
fact, different embedded monopole content. 

We also study the relative ratio R of the monopole density in the IR clusters pjfj 3 compared 
to the total monopole density , 


R 


„latt 
RIR 
n latt 
r total 


(44) 


In order to extrapolate this ratio to the continuum limit we use the linear formula: 


fl“( a) = /,■„('I - Ka). 


(45) 


Here Rq and K are the fitting parameters. An example of the extrapolation and the extrap¬ 
olated values of R are shown in Figures |H1 (a) and (b), respectively. Note that according to 
Eq. the formula for extrapolation EH should contain 0(a 2 ) corrections, which, however, 
can not be traced out due to limited accuracy of our data. 


D. Discussion on scaling properties 

It is interesting to speculate about the nature of the observed scaling behavior of the 
embedded monopole densities (EH), EH, EH, EH, EH- In the zero temperature case 
the naive physical density, p phys (a, A) = a _3 p latt (a, A), of the scalar, axial and invariant 
embedded monopoles diverges in the continuum limit as 

pglK A) ~ a -1 , pf ys (a, A) ~ a ~ 2 , A ^ 0, (46) 

pg7(a,0) ~ a" 1 , pj hys (a, 0) ~ gU 3 , A = 0. (47) 

Let us suppose for a moment that the general quantity EH is a density of objects of an 
unknown dimension, and that the individual objects are not strongly correlated in the lattice 
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FIG. 8: (Color online) (a) Example of extrapolation of the infrared-to-total ratios (1441) correspond¬ 
ing to the scalar, axial, and chirally invariant embedded monopoles for the eigenvalue A = 235 MeV. 
The fit is done by the linear formula (1451) . (b) The infrared-to-total ratios of the scalar, axial and 
invariant monopole densities extrapolated by Eq. (03) to the continuum limit vs. A. 


ensembles. Then, on general grounds, the terms in the formula m can be interpreted as 
follows. 

If the object is of a pure lattice origin (a lattice artifact), then its lattice density should 
not change with the variation of the physical scale a. Thus, if the first term C is non-zero 
in the continuum limit a —> 0, then density p corresponds to a purely lattice object. The 
physics of these objects is determined by the ultraviolet cutoff Auv ~ a” 1 only. 

Now suppose, that parameter the C vanishes and the leading behavior of the lattice 
density is pi att = va + ... as a —> 0, where the coefficient v is of the order of the physical 
QCD scale v ~ Aqcd- Then the world-manifolds of the objects are the three-dimensional 
volumes distributed in the four-dimensional space time with the physical density v. The 
corresponding object is a membrane. 

The leading scaling behavior in the form pi att = sa 2 + ... corresponds to string-like objects, 
density of which is given by the quantity s ~ Aq CD . Finally, if one studies pointlike objects 
(he., monopoles), which are not strongly correlated, then the scaling of their density should 
be pi a tt = pa 3 + ..., and the physical density of the objects should be p ~ Aq CD . 

However, these simple considerations become incorrect if the objects are (strongly) in¬ 
teracting with each other. As an illustrative example, it may be useful to consider currents 
of Abelian monopoles in an Abelian projection of pure SU(2) Yang-Mills theory. A general 
configuration of the gauge fields typically contain 2f| 0] two components of the monopole 
clusters, one of them is infrared and the other one is ultraviolet. The physical density of the 
infrared monopole currents is finite in the continuum limit a —> 0, which means that the for 
these monopoles the coefficients C, v and s in Eq. ED are zero. On the other hand, the 
ultraviolet component of the Abelian monopole density diverges as a~ l in physical units. 
One can understand this scaling as a consequence of a strong correlation between segments 
of the monopole loops at the scale of the lattice spacing a because a typical UV monopole 
cluster is, in fact, a loop of the length of a few lattice spacings. One can equivalently say 
that the monopole clusters are short-ranged dipoles. The nature of this strong correlation 
is of a purely lattice origin as the recent data shows 27]. Indeed, it was found in Ref. 27] 
that the density of the UV monopoles strongly depends on the UV-properties of the gluon 
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action. The density of the IR monopole clusters are also sensitive to the lattice details of 
the gluon action, since the artificial UV monopoles may randomly connect to the physical 
IR clusters and be counted by a lattice algorithm as a part of the physical IR cluster. 

Thus the a -1 and a~ 2 scaling of the physical densities of the scalar/axial and, respectively, 
invariant embedded monopoles may be a result of the lattice procedure(s) which may be 
sensitive to the UV-scale. On the other hand one can not exclude a possibility that the 
embedded QCD monopoles may be strongly correlated with objects which have surface-like 
and 3D volume-like world trajectories. This property is supported by the observation j2^| 
that the low-lying fermion modes show unusual localization properties being sensitive both to 
the physical scale Aq C d and to the ultraviolet cut-off, a -1 . If this suggestion is correct, then 
the scaling of the “slave” monopoles may manifest the scaling of the “master” objects. In 
this case the scalar/axial and invariant monopoles should be correlated with (or, as one can 
also say, “lie on”) strings and membranes, respectively. In this paper we are not performing 
a detailed scaling analysis of the monopole clusters concentrating on simplest properties 
only. A review of the lattice data on many-dimensional vacuum objects in four-dimensional 
Yang-Mills theory can be found in Ref. |28j. 

Finally, let us note that the embedded monopoles are likely not lattice artifacts because 
their densities scale as non-zero powers of the lattice spacing a. Still, the effect of the lattice 
artifacts on the density may be noticeable since the gauge fields are not improved at all 
contrary to the Dirac eigenmodes which are greatly improved by using the overlap fermions. 


IV. EMBEDDED MONOPOLES AND FERMION SPECTRAL DENSITY 

One of the most essential characteristics of the fermion modes in the gauge theory is the 
fermion spectral density pp which is formally defined as the expectation value 

PfW = \ jf£>(A-A)>, (48) 

A 

where the sum goes over all Dirac eigenvalues A = A(A) corresponding to gauge fields config¬ 
urations A which enter the partition function. The low-lying part of the fermion spectrum 
is important for the chiral symmetry breaking due to the Banks-Casher formula [30j, 

(W*) = ~ lim7rp F (A), (49) 

A—>0 

which relates the chiral condensate (ipip) with the spectral density. 

Examples of the spectral density pf as functions of A are shown in Fig. El for the con¬ 
finement (T — 0, (3 — 2.3493) and for the deconfinement (T = 1.15T C , (3 = 2.35) phases. 
In the zero temperature (confinement) case the spectrum is a gradually increasing function 
of the eigenvalue A, and the low-energy part of the spectrum has a finite A —> 0 limit. In 
the deconfinement phase the low-energy part of the spectrum is suppressed compared to the 
confinement phase but is, however, non-zero. This feature as well as the peak of the spectral 
density at A = 0 may be an artifact related to an effect of the finite volume. We expect that 
in the limit of an infinite volume the spectrum above the deconfinement temperature should 
vanish below some critical value A C (T). This property could imply vanishing of the chiral 
condensate © in the deconfinement phase, > T c ) = 0. In our case the critical 

value is A C (T = 1.15 T c ) « 400 MeV. 
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FIG. 9: (Color online) Spectral fermion density in the confinement (T = 0, (3 = 2.3493) and in the 
deconfinement (T = 1.15T C , /3 = 2.35) phases. 

The embedded QCD monopoles are suggested |3] to be agents the chiral symmetry 
restoration, since in the cores of these monopoles the chiral invariance should be unbro¬ 
ken. According to the proposed scenario, one can expect that at low Dirac eigenvalues - 
which are relevant to the chiral symmetry breaking due to the Banks-Casher relation © 
- the density of the embedded monopoles should be high in the chirally invariant (high 
temperature) phase and the density should be relatively low in the chirally broken (low 
temperature) phase. Thus suggestion implies, in turn, that the density of the embedded 
monopoles should be anti-correlated with the fermion spectral function: the lower value of 
the spectral function the higher monopole density is expected to be. In particular, in the 
high temperature phase the vanishing spectral function at A < A c implies the high density 
of the quark monopoles for A < A c , and vice versa. 

The anti-correlation of the quark monopole density and the fermion spectral function is 
indeed observed in deconfinement phase as it is shown in Figure ^d). Indeed, as one can 
see from comparison of Figure^d) and the T = 1.15T C spectral function shown in FigureEI 
at high A the fermion spectral function is high corresponding to low embedded monopole 
density, while at low A the fermion spectral density is suppressed in accordance with the 
observed large valued of the monopole density. 

In the confinement case the qualitative relation between the spectral density and the 
embedded monopole density is true as well according to Figures IHa),(b),(c) and Figure El 
the fermion spectrum is an increasing function of the Dirac eigenvalue A while the monopole 
densities are generally decreasing function of A. 


V. EXCESS OF GLUON ACTION ON QUARK MONOPOLES 

If the embedded QCD monopoles are physical objects then we would expect that these 
objects are locally correlated with the action density and, presumably, with the topological 
charge. Following Refs. mm we calculate numerically the excess of the SU(2) gauge 
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action at the position of the embedded monopole currents j X4n 

f {\jx,/l\ ScX'/j) — 6(\jx,fl\) (Sp) /r n \ 

Js = -TP—rc-i j (Y U J 

where S c is the sum over the elementary plaquette actions, 

Sc=J2 Sp ’ Sp = l-^TvU P , 

Pedc 

belonging to the six faces P of the cubes c = c x ^ with non-zero monopole charge, j Xtfl ^ 0. 
Here Up is the SU(2) plaquette constructed from the lattice links U Xfl in the standard way 
Up x ^ = U xlx U x+ jx^U' xJri>4J Ul v . The second term in Eq. (I5U1) subtracts the vacuum average of 
the action from the value of the gluon action at the monopole. 

Equation © can be understood as the (average) excess of the Yang-Mills action calcu¬ 
lated at the (average) distance r = a/2 from the center of the monopole. In fact, in any 
given lattice configuration the position of the monopole center can not be determined ex¬ 
actly within the lattice the cube possessing a nonzero monopole charge. However on average 
the monopole center is located as the cube center which resides at the distance a/2 from 
any face (plaquette) of the cube. It is worth noticing that Eq. (I5U1) defines the excess of the 
chromomagnetic part of the action since by construction the action around the monopole 
is calculated on the plaquettes P perpendicular to the corresponding link l = { x , /a} of the 
monopole trajectory. 

In the naive continuum limit the elementary plaquette, say P x ,i 2 , is expanded in powers 
of the lattice spacing as follows 


Sp,M =a“^ + Op) 


(51) 


where g is the (bare) coupling and U XI1 = e ia 9 A U x +p / 2 ). Thus the excess of the action (1301) 
can be written as (we omit 0(a 6 ) corrections starting from here) 


where 


7T 


f s (a) = a 


a/2), 


(52) 


B, m (a/2) = (a,(B') 2 ), 


r=a/2 


- (o, (B c ) 2 >, 


(53) 


is the excess of the chromomagnetic condensate at the distance r = a/2 from the quark 
monopole and a s = g 2 j[\.T\). The chromomagnetic field at the segment of the monopole 
current j y is defined as B/(j u ) = ~e flua gF/ l g/2. This definition reduces to the standard 
one for a static (v = 4) monopole: B/ = e i ^ k F ( / k /2 with i,j,k = 1,2,3. In the Euclidean 
space-time at zero temperature the chromomagnetic and the standard gluon condensates 
are related as ( a s (B c ) 2 ) = (a s (F/ u ) 2 )/2. 

Before proceeding with analysis of the numerical data it would be instructing to discuss 
the expected behavior of the chromomagnetic fields inside the embedded monopoles at least 
on a qualitative level. In Ref. Jl] the embedded QCD monopole is associated with the Narnbu 
monopole in the Electroweak model. The Narnbu monopole is essentially the ’t Hooft- 
Polyakov in na (HP) monopole configuration embedded into the EW model. Therefore 
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one naively can expect that the behavior of the chromomagnetic fields inside the embedded 
monopole in QCD is qualitatively similar to that of the HP monopole in the Georgi-Glashow 
model. 

As an illustrative example let us consider the Bogomohny-Prassad-Sommerfeld (BPS) 
limit [HI of the Georgi-Glashow model, 


A 


2\2 


£gc = - a (F;„) 2 + - (D?<ty + - ((4“) 2 - r, 2 ) 


(54) 


This model describes the dynamics of the SU{2) gauge held A“ interacting with the triplet 
(adjoint) Higgs held <h a , a = 1, 2, 3. The adjoint covariant derivative is given in Eq. (J7J). The 
scalar coupling A describes self-interaction of the Higgs held. The condensate of the Higgs 
held is |($)| = rj. The masses of the gauge and Higgs helds in the Georgi-Glashow model 
are, respectively, itia — gv and m$ = V2Xg. 

The BPS limit is defined by the condition A = 0, which sets the mass of the Higgs particle 
to zero, = 0. Due to the absence of the quartic Higgs self-interaction the classical static 
’t Hooft-Polyakov solution can be found explicitly 


34 : 


= — 2 H{gg r), A? = e ai ~[l - K( V g r )], = 0, (55) 

gr 2 gr 2 


where 


m = 


sinh £ ’ 


ff«) = {coth5-l, 


(56) 


The chromoelectric held of the HP monopole is zero, Fj* = 0, and the chromomagnetic held 
B t = e ijkFj k /2, is: 


m = 


g f 


K 2 - HK) + —- HK , 
/ g r 2 


(57) 


The corresponding “condensate” of the chromomagnetic held tends to a finite value in the 
monopole center, r —> 0: 


BhpM . g[B 2 W ] 2 = ^ 


B„ P (r) = B 


(0) 

HP 


1 + O ( (rjgr ) 


1 — K 2 (r]gr)J + 2H 2 {r]gr) K 2 (rjgr ) 

„4 „4 

(o) _V 9 for r ^ 0i 


K' u > — 
°HP — 
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(58) 

(59) 


Let us take for a moment this illustrative example of the hedgehog configuration seri¬ 
ously. The QCD counterpart of the held d>“ is an octet quark-antiquark composite held £ a . 
We do not expect the presence of the octet condensates (£ a ) in vacuum of the Yang-Mills 
theory because such a condensate must inevitably break the color symmetry (see, however, 
a discussion in Ref. [Hj]). On the other hand one may expect [jj] that the non-perturbative 
color-invariant four-quark condensates ([36] of the form ((£ a ) 2 ) should stabilize the hedgehog¬ 
like configurations made of the composite “Higgs” held £ a in the confinement phase. 

One can expect that the value of the condensate rj in Eqs. (EH), (EH) should be of the 
order of a typical dimensional quantity describing the chiral condensate, g ~ I W)| :1/3 
0.2 ... 0.3 GeV, which, in turn, is of the order of the QCD scale parameter Aqcd- One can 
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think of g as of the condensate outside the core. The gauge coupling g can be associated 
with the QCD running coupling. Then Eq. (EH) predicts that the chromomagnetic field 
inside the embedded monopole should be “soft”, 

Bmon (r) ~ 5 ,4 ( r A qcd) • Aqcd • 


In particular, this relation implies the absence of the hard ultraviolet divergences of the 
energy density inside the monopole cores contrary to the case of the Abelian (Dirac) 
monopoles jdj. E/]. 

We perform the fit of the numerical data for the correlation function dsn by 


(•fit,HP 
JS 



+ 2 K 2 {x)H' 2 (x) , 


x = g(Anp a/2) • rja/2 , 


(60) 


which can be obtained from Eqs. '& dSHD by identihcation r = a/2. The function K is 
given in Eq. (ESI), and instead of the Georgi-Glashow coupling g we take the one-loop running 
coupling constant of the SU(2) Yang-Mills theory, 


g 2 (A a) 


11 

12t r 2 



(61) 


The htting parameters are the HP scale parameter Ahp and the “condensate” parameter g. 




FIG. 10: (Color online) The excess of the chromomagnetic action (15011 at (a) the scalar and (b) 
the chirally invariant monopoles at A = 0 and A = 235 MeV. The solid and the dotted lines refer 
to the fits Eq. m and m, respectively. 

Examples of the correlation function fs at A = 0 and A = 235 MeV are shown in Fig¬ 
ures m a) and (b) for scalar and chirally invariant monopoles, respectively. The excess of 
the chromomagnetic energy for the scalar and axial monopoles coincide with each other 
within error bars. The examples of the fits of the excess energy by the function dUD are 
shown by the solid lines. One can see that the fitting formula Egp, dnnp , ed - which is 
resembling the’t Hooft-Polyakov monopole configuration in the Bogomol’ny limit ESI), ESI), 
ESI) - emulates our numerical data relatively well. The fitting is better near the continuum 
limit, a —> 0. However at relatively large distances, a > 0.13 fm, the data and the best fit 
function show some noticeable difference which makes y 2 /d.o./ ~ 3... 5 for these fits. 

The best fit parameters Ahp and g obtained in our fits by the function ESD come with rel¬ 
atively large errors. For example, for scalar monopoles at A = 0 we have Ahp = 17(26) MeV 
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and g = 136(14) MeV, while at A = 235 MeV we get A H p = 28(31) MeV and g = 98(12) MeV. 
The corresponding numbers for the invariant monopoles are: in the A = 0 case we obtain 
A H p = 30(36) MeV and g = 72(17) MeV, while for A = 235 MeV we get A H p = 24(20) MeV 
and g = 75(10) MeV. Thus the values of the parameter A H p can not be defined well due to 
quite weak dependence of the logarithm function © on the value of its argument. More¬ 
over, the value of App is very small what makes g ~ 1 in all our fits. On the other hand 
the fit quantitatively confirms that the values for the “condensate” 77 is of the order of the 
QCD scale, 77 ~ Aqcd- The effective size of the “HP monopole” core, 1/rriA = 1 /(gg), can 
be estimated to be of the order of 1 fm for all values of A. Since this value is unrealistically 
large we conclude that the fact that the HP fitting formula © works relatively well is 
only a manifestation of the “softness” of the gluonic action inside the core of the monopole. 
Quantum corrections to the HP monopole fields may be important. 

In order to get prescription-independent result on the scaling of the average action excess 
we fit the available data by the power-like fitting function 

/J*-p° wo » = (a/h ) a+s , (62) 

where the scale h and the “anomalous” exponent 5 are the fitting parameters. If the action 
density is independent on the distance to the monopole center (or, in other words, if the core 
of the embedded monopoles is structure-less) then we would naively expect the vanishing 
anomalous exponent, 5 = 0. If 5 7 ^ 0 then one can expect some structure of the monopole 
core. 

The example of the fits © at A = 0 and A = 235 MeV are shown in Figures HUT al and 
(b) by the dotted lines. As one can see from these Figures, the HP monopole fit (1HU1) and 
the simple power fit © are practically indistinguishable from each other. Note that both 
fits are two-parametric ones. 




FIG. 11: (Color online) The best fit parameters (a) h and (b) 5 of the fit function (16211 vs. the 
Dirac eigenvalue A for the scalar, axial and chirally invariant embedded monopoles. 


The best fit parameters for the power function © are shown in Figures ©a) and (b) 
as functions of A for all studied types of the quark monopoles. We find that both h and 5 
parameters are almost independent on the Dirac eigenvalue A. Moreover, these parameters 
for different types of the quark monopole are quite close to each other. This result suggests 
that the gluonic structure of the embedded monopoles seems to be independent on the value 
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of A. The typical values of the fitting parameters are concentrating around central values 
hg A ~ 1.6 fm, hi rs 2 fm and Ss,a,i ~ 0.7 with, however, relatively large error bars. 
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FIG. 12: (Color online) The local excess (1501) of the Yang-Mills action at the position of the (a) 
scalar and (b) invariant monopoles, scaled by the power function, /g(a) a -( ' 4+< b, with the anomalous 
scaling exponent <5 = 1/2. The three values of the lattice coupling /3 are shown. 


In order to get an impression on the quality of the power-like scaling © we plot in 
Figures H21 (a) and (b) the scaling function fs multiplied by the factor a~^ +s ^ with 5 = 1/2 
for the three different values of the lattice coupling constant [3. Figures H2K a) and (b) clearly 
show that the quantity / 5 (a) a -4,5 is independent of the lattice spacing a for the embedded 
monopoles of all three types. Since the fits of by the power © and the HP-inspired © 
functions are practically indistinguishable at our data sets (as one can see from examples 
plotted in Figures HUT a) and (b)), the scaling of our data with the HP-type energy excess (l60l) 
should be as remarkable as it is plotted in Figures m for the power-profile excess ©. 

It is interesting to point out that despite noticeably different scaling (1331) of the 
scalar/axial and chirally invariant monopole densities towards continuum limit, the excess 
of the action density at the positions of these monopoles scales essentially in the same way. 

Since the monopole cores are “soft” we do not expect a fine-tuning Y7] between the 
energy and entropy of these monopoles (at least, in the studied case of the pure SU(2) 
Yang-Mills theory). On the other hand, the embedded monopoles does not scale as particle 
like objects, therefore the investigation of the energy-entropy balance in the case of the 
embedded monopoles in the quenched case may be a complicated issue. 

A cautionary remark here is that the monopoles are detected with the help of the oper¬ 
ators m which implicitly depend on the lattice spacing a. If the monopole has a core-like 
structure of the size of the lattice spacing then the ability of the numerical procedure to 
“detect” a lattice monopole within the lattice cube should be very sensitive to the size of 
the “detector” (be. to the lattice spacing). I 11 order to get physically reliable results on the 
monopole density and the monopole correlations one should probably study the scaling of 
the extended (blocked) monopoles [38]. 

Summarizing, both scalar, axial and invariant embedded monopoles are locally correlated 
with the magnetic part of the gluonic action. The excess of the action at fixed lattice spacing 
a (or, equivalently, at fixed (3) is a slowly increasing function of the Dirac eigenvalue A. The 
average action excess on the scalar and axial monopoles coincide with each other and is 
approximately two times higher than the excess of the action on the invariant monopoles. 
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The positive value of / 5 (a) indicates that the action density near the monopole is increased 
compared to the average density. Thus the embedded monopole has a chromomagnetic core. 
On the other hand, the positive value of the anomalous scaling exponent 6 indicates that the 
excess of the chromomagnetic action decreases as one approaches the center of the monopole 
core. Concretely, the chromomagnetic condensate vanishes in the center of the embedded 
monopoles as 


£mon(r)~ r 5 , 5 ~ 0.5 ... 1. (63) 

Thus, the embedded monopoles possess “chromomagnetically” empty cores. 


VI. CONCLUSIONS 

We study basic properties of the embedded QCD monopoles in the pure SU(2) Yang- 
Mills theory. The monopole trajectories are found with the help of the low-lying eigenmodes 
of the overlap Dirac operator. These modes are then treated as c-valued quark fields, the 
behavior of which emulates the basic chiral properties of the QCD vacuum. 

The embedded monopoles are explicitly gauge-invariant, and the magnetic charge of the 
monopole is quantized and conserved. Basically, the embedded QCD monopoles are the 
gauge-invariant hedgehogs in the quark-antiquark condensates (therefore we also call these 
monopoles as “quark monopoles”). 

We give the lattice definitions of the embedded monopoles of scalar, axial and chirally 
invariant types. We find that the scaling of the scalar and axial monopole densities towards 
the continuum limit is the same as the scaling of the string-like objects. The scaling of 
the chirally invariant monopoles corresponds to the one of the membrane-like objects. This 
result may indicate that the monopole trajectories are correlated with higher-dimensional 
(string-like and membrane-like) objects in the SU( 2) Yang-Mills theory. The “scalar/axial 
string”, for example, may be a border of the “chirally invariant membrane”. We also observe 
a difference in the scaling properties of the monopoles corresponding to the non-zero and to 
the zero Dirac eigenvalues. 

The embedded QCD monopoles were suggested [ 3 ] to be related to the restoration of the 
chiral symmetry in the high-temperature phase since their cores should contain the chirally 
symmetric vacuum. Our numerical study supports this suggestion since the monopole den¬ 
sity is anti-correlated with the density of the Dirac eigenmodes. In particular, at low Dirac 
eigenvalues - which are relevant to the chiral symmetry breaking due to the Banks-Casher 
relation - the density of the embedded monopoles is high in the chirally invariant (high 
temperature) phase and is relatively low in the chirally broken (low temperature) phase. 

We find that the embedded monopoles have gluonic cores, which are more pronounced for 
the chirally invariant monopoles compared to the scalar/axial monopoles. On average, the 
chromomagnetic energy near the monopole trajectories is higher compared to the chromo¬ 
magnetic energy far from the monopole core. However, our scaling analysis suggests that at 
the very center of the embedded QCD monopole, the excess of the chromomagnetic energy 
reduces back to the vacuum expectation value. Therefore a typical monopole core is a bump 
in the chromomagnetic energy which takes its maximum value at a certain finite distance 
from the monopole center. Outside this bump - towards the monopole center and/or far 
from the monopole core - the energy density diminishes to its vacuum expectation value. 
This structure is similar to the structure of the ‘t Hooft-Polyakov monopoles if one attributes 
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to the asymptotic freedom the suppression of the chromomagnetic gluon condensate in the 
monopole center. 

Finally, we would like to remark that one can not exclude a possibility that the properties 
of the embedded monopoles in the full QCD may drastically be different from the quenched 
case. On the other hand the quenched theory mimics chiral instability of the full QCD 
to develop a chiral condensate at low temperatures. Therefore our results support the 
suggestion that the quark monopoles are tightly related to the chiral symmetry restoration 
also in the case of the real QCD. 


Note added 

When numerical calculations were completed we became aware that a similar procedure 
has been applied earlier in SU(3) lattice gauge theory (39|. 
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APPENDIX A: DETAILS OF NUMERICAL SIMULATIONS 

We simulate numerically the SU(2) lattice gauge theory with the standard Wilson action, 
Sp = (3(1 — (l/2)Tr Up), where (3 is the SU(2) gauge coupling and Up is the SU(2) plaquette 
variable constructed from the lattice link fields U(x,n). We used various values of f3 at 
different lattice sizes to check the scaling of the numerically calculated observables towards 
the continuum limit. The parameters of the numerical simulations are given in Table' Q 


T = 0 

P 

a[fm] 

L s 

Lt 

N C onf 

P 

a[fm] 

L s 

L t 

N C onf 

P 

a[fm] 

L s 

L t 

Nconf 

2.3493 

0.1397(15) 

10 

10 

300 

2.3772 

0.1284(15) 

10 

14 

90 

2.4071 

0.1164(15) 

12 

12 

180 

2.4180 

0.1120(15) 

12 

14 

150 

2.4500 

0.0996(22) 

14 

14 

200 

2.5000 

0.0854(4) 

16 

18 

200 

T = 1.15 T c 


2.3500 

0.1394(8) 

16 

4 

200 



TABLE I: (Color online) Parameters of the simulations. 


The first 6 points in Table [0 correspond to the zero temperature (confinement) phase. 
The lattice geometries and values of the lattice coupling (3 are tuned in order to keep the 
lattice volume constant, V = 3.8 fm 4 . The point with /3 = 3.5 has a little bigger volume, 
V = 3.92 fm 4 . 
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In order to have an impression about the behavior of the quark monopoles in the high 
temperature phase we study one point at asymmetric lattice 16 3 x 4. At these lattices the 
system is just above the finite temperature critical point with T = 1.15 T c . 

In order to define the quark monopoles one may use eigenmodes of the Dirac operator 
in the background of the gauge held. In our lattice simulations we use the overlap fermions 
which possess an exact chiral symmetry, enjoy automatic 0(a) improvement, and are not 
defaced by exceptional configurations [4(j|. The overlap Dirac operator is 

D = ^(i + D w /Jd w dI) =£(1 + 75 sign (if)), H = lb D w , (Al) 

where D w is the Wilson Dirac operator with negative mass term and H is hermitian Wil¬ 
son Dirac operator. The value of p parameter is equal to 1.4. We have used the minimax 
polynomial approximation to compute the sign function. In order to improve the accuracy 
and performance about one hundred lowest eigenmodes of H were projected out. The eigen¬ 
values of D , which lies on the circle in the complex plain, were stereographically projected 
onto the imaginary axis in order to relate them with continuous eigenvalues of the Dirac 
operator. 
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